Symmetry, Integrability and Geometry: Methods and Applications 



Vol. 2 (2006), Paper 065, 15 pages 



On the Linearization of Second-Order Differential 
and Difference Equations 

Vladimir DORODNITSYN 

Keldysh Institute of Applied Mathematics of Russian Academy of Science, 
4 Miusskaya Sq., Moscow, 125047 Russia 
E-mail: dorod@spp.keldysh.ru 



Received November 28, 2005, in final form July 13, 2006; Published online August 16, 2006 
| Original article is available at | |http : //www . emis . de/ j ournals/ SIGMA/2006/Paper065/| 

^3 ■ Abstract. This article complements recent results of the papers [J. Math. Phys. 41 (2000), 

<^ | 480; 45 (2004), 336] on the symmetry classification of second-order ordinary difference 

equations and meshes, as well as the Lagrangian formalism and Noether-type integration 
technique. It turned out that there exist nonlinear superposition principles for solutions of 
special second-order ordinary difference equations which possess Lie group symmetries. This 
superposition springs from the linearization of second-order ordinary difference equations 
\yy by means of non-point transformations which act simultaneously on equations and meshes. 

These transformations become some sort of contact transformations in the continuous limit. 
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O ' 1 Introduction 

oo " 

The recent article was devoted to symmetry classification of second-order ordinary difference 
■ equations. This was based on the paper by S. Lie [Sj, in which he has provided a symmetry 

classification of second-order differential equations (ODEs). A special class of the specified in- 
variant second-order difference equations converges to Lie's invariant second-order ordinary 
differential equations in the continuous limit. The Lagrangian formalism for Lie's list of those 
. second-order ordinary difference equations associated with special invariant meshes (lattices) 

was developed in recent paper [2]. From the Noether theorem [2] it is known that invariance of a 
variational functional with respect to an r-parameter group leads to existence of r conservation 
laws for the corresponding Euler's equation. It was shown in [2] (see also H>] for details) 
that there exists a similar (although more complicated) construction for difference models. This 
discrete analog of the Noether theorem can be applied to analytic integration of difference equa- 
tions. Namely, two-dimensional symmetries of the Lagrangian for a second-order difference 
equation provide two first difference integrals of the two-point type. This yields a general 
solution of the difference system by purely algebraic manipulations. The last point is impor- 
tant for difference equations since there are no analytic integration techniques for constructing 
solutions. 

The problem of linearization of a given equation (or a system of equations) is of permanent 
interest in mathematics. For example, in the paper [Jj necessary and sufficient conditions, 
under which a given nonlinear system of PDEs with n > 2 independent variables and m > 1 
dependent variables can be transformed to a linear system of PDEs, were developed. It was 
shown that in all the cases of such mappings, an infinite-parameter Lie group of point or contact 
transformations admitted by the original nonlinear system is needed. 
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For ODEs the situation is different. S. Lie proved that symmetry algebras admitted by 
second-order ODEs are finite, and the maximum symmetry algebra is 8-dimensional. This max- 
imum symmetry algebra is admitted by a linear equation. S.Lie has developed necessary and 
sufficient conditions under which a given nonlinear second-order ODE can be transformed into 
a linear one jSJ (the discussion of this question can be found, for example, in p. 38]. One 
of these (equivalent to each other) conditions is existence of two linearly connected symmetry 
operators. In this case there exists a change of variables that transforms a given nonlinear 
equation into a linear one. Generalization of this idea to ordinary difference equations was de- 
veloped in JU] where, for some class of ordinary difference equations on fixed uniform lattice, 
point transformations which linearize corresponding nonlinear difference equations have been 
found. This class of ordinary difference equations was restricted by those equations that pos- 
sess point autonomous symmetries, which do not change the independent variable. Difference 
equations are assumed to be defined on a fixed regular lattice with unit spacing. This ap- 
proach was generalized in jllj to non-autonomous symmetries, while the independent variable 
was still fixed. The discrete analog of Lie's theorem jSj, providing the necessary and sufficient 
conditions, under which a given nonlinear second-order difference equation and a mesh can be 
transformed into a linear one, was developed in Theorem 1]. It is important to notice that 
in all the cases the corresponding point transformations change equations within one class of 
equivalent difference equations, which are invariant under the actions of isomorphic symmetry 
groups. Point transformations of such type cannot connect nonlinear equations from Lie's list 
with a linear one. Being based on listing of groups, the symmetry classification of second-order 
differential [2j and difference equations singled out classes of equations, which cannot be 
connected by means of point transformations (of independent and dependent variables) and, in 
particular, cannot be linearized by a point transformation in principle. In the classification [2j 
S. Lie singled second-order ODEs up to arbitrary point transformations. For example, a linear 
equation is the simplest representative equation from the class of equivalent equations which 
possess an 8-dimensional symmetry group. 

Thus, in this article we are dealing with "linearization of not-linearizable equations" in the 
above mentioned sense. We establish non-point transformations connecting some nonlinear 
invariant second-order differential equations from Lie's list with linear equations. Such trans- 
formations have discrete counterparts, and so the invariant difference schemes can be linearized 
together with difference meshes. 

We will use the following notation. Let x be the independent variable and y the depen- 
dent variable. To consider a second-order difference equation and a lattice in x-direction we 
need 3-points of a difference mesh (lattice): (x,X-,x + ) and corresponding values of the de- 
pendent variable in three neighboring points: (y, y_, y+), where "+" marks a point shifted to 
the right, and "— " marks a point shifted to the left. We call this set of points (a subspace) 
(x,x^,x +1 y,y^,y + ) a difference stencil. It is supposed that the stencil and the corresponding 
difference equation can be shifted to any other point of the lattice. 

The discrete model for second order ODE can be presented in terms of two difference equa- 
tions: 



We will also use the following notation for the right and for the left difference derivatives of 
the first order: 






y+ - y 



y-y- 
h- 



y x = 



ys = 



where h+, h- are the space steps of the mesh in the x-direction. 
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To distinguish between continuous derivatives and difference ones we denote the first by 
y',y", ■ ■ ■■ The first equation of the system represents a second-order difference equation. 
In the continuous limit it should become a second-order ordinary differential equation on some 
point of a finite-dimensional space of continuous variables. The second equation yields a lattice 
(mesh) on which the first equation is considered. Notice that the mesh equation does not need 
to be a second-order differential equation in the continuous limit. For example, one can consider 
a uniform lattice equation h + = h— (where h + = x+ — x, h- = x — X-) which "disappears" in 
the continuous limit. 

The group transformations considered in the present approach are of the same type as for 
ODEs. They are generated by a Lie algebra of vector fields of the form 

d d 

x = ^ x ' y) dx~ + r]{ - x ' y) dy~- 

The corresponding transformations are purely point ones, since the coefficients £ and n 
depend on (x,y) only, and do not depend on the shifted points (x + ,y + ) or (x_,y_). 

In the classification of second-order ODEs [3] S. Lie identified a linear equation, which possess 
8-dimensional symmetry group; four simplest representative nonlinear equations which possess 
a 3-dimensional symmetry group; 2 classes of ODEs with a 2-dimensional symmetry group and 
one class of ODEs with a one-dimensional symmetry group. 

In the following sections we consider three examples of ODEs that possess three symmetries. 
In parallel we will consider the corresponding difference models as developed in Also we 
consider examples from the Lie classes of nonlinear ODEs with a two-dimensional Lie algebra. 
The chosen equations are important from the physical point of view. 

We will write down explicitly the linearizing tangent transformations for underlining second- 
order nonlinear ODEs. We will also write down explicitly the linearizing two-points transfor- 
mations for corresponding difference equations and meshes. Notice that these non-point trans- 
formations change the admitted symmetries of the underlying equations as the corresponding 
linear equations possess the 8-dimensional symmetry. In fact, these transformations are not 
explicitly connected with admitted symmetries of ODEs. This aspect will be discussed later. 

Notice that all difference and differential equations considered here are integrable (at least in 
quadratures) as far as they possess a symmetry algebra of dimension not less than 2. Therefore, 
the computational aspects of the invariant difference schemes (as well as stability, convergence 
of developed schemes, etc.) will not discussed here. 

In the present paper we are following the set of invariant ODEs for which an appropriate 
invariant difference equations and meshes were developed in PQI2J- 

It seems that the main result presented below is the explicit non-point transformations, which 
linearize second-order ODEs from the Lie's list and appropriate invariant difference equations, 
which possess symmetries of the dimension less than 8. A discussion of the origin of such 
transformations will be presented at the end of Example 1 and in Concluding remarks. 

2 Example 1: ODE with a 3-dimensional symmetry group 

We start with the example of a second-order nonlinear ODE from the Lie's list of equations, 
which possesses three symmetries. It has two first integrals and, therefore, can be integrated by 
means of algebraic manipulations with difference integrals. 
The group given by the operators 
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corresponds to the invariant differential equation from the Lie's list, 

y" = y~ 3 - (3) 

This equation can be obtained from a variational functional with the Lagrangian function 
L = y' 2 — 4^ which admits all three operators © as variational symmetries. 
The Noether theorem 0] yields the following first integrals: 

Ji = y' 2 + \ = a , h = 4 - (v - y' x )y' = B °- ( 4 ) 

yZ yZ 

By means of the first integrals we write the general solution y(x) as 

A y 2 = (A x - B ) 2 + 1. (5) 

Now we reproduce the procedure of the linearization for the equation (j2J), followed by ex- 
planatory comments. 

It turns out that equation (J3J) can be linearized by the non-point transformation 

A , 1 A , „„^„„^.„.n „.„./ ..t„\ 1 „./ ../ 



t = x [y' + -g + arctan(yy') - yy', y{x) = — — , y' = -u', (6) 

which transforms equation © into the following linear equation 

u" + u = 0, (7) 

which possesses a linear superposition for its solutions. 

Thus, the developed formulas implicitly reproduce the nonlinear superposition principle for 
equation (|5]). The transformation © changes the first integrals (j3J) into the following ones with 
the same constants Aq,Bq: 



v! + u = Aq, t + arctan I — I = Bq. 

The latter yields the general solution of the linear equation 

u(t) = \/~Aq cos(Bq —t) = A sin t + B cos t . 

Notice that (jUJ connects the general solution © with (10), which provides an explicit proce- 
dure for the superposition of solutions (JSJ). Indeed, let us take two solutions among the family (J5J), 
which correspond to the two sets of constants: y\ = Y\{x, Ai, B\) and y% = Y2(x, A2,B2)- Then 
the transformation (jHJ) produces two solutions of the linear equation (J7J), 

m(t) = y/A\cos(B 1 - t), u 2 {t) = v^cos^ - t). 

The sum of these solutions 

u(t) =ui+u 2 = \/~A*cos(B* - t), 

^AisinBi + vC42sini?2 



A* = A\ + A 2 + 2y/A 1 A 2 cos(B 1 - B 2 ), B* = arctan 



' A\ cos B\ + V A 2 cos B 2 

yields a new solution of equation (j3J) with two constants A*, B*. The developed formulas 
reproduce the nonlinear superposition principle within the family of solutions ©• For example, 



the solutions y\ = V x 2 + 1 and y 2 = J Ax 2 + j produce the new one ?/3 = J 7x 2 + j . 
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Now we consider the discrete case. In the paper [2] a difference scheme which preserves 
symmetries and the Lagrangian structure of the initial equation © was developed. 
The starting point was an entire set of finite-difference invariants 

t h + t h - T y 2 y- (v+-y v-v- 



yy+ yy- h_ \ h + h- 

in a subspace of difference variables (difference stencil) x, h + , h_, y, y + , y_. 

Any difference equation which can be represented by means of the above invariants will be 
invariant. Thus, in such a way one can generate an invariant mesh. We will use the following 
relation between difference invariants 

(8) 



yy+ yy- 

as the invariant mesh (see j2j). 

It is easy to see that the mesh equation has a first integral 

e, e = const. 



yy+ 

As an invariant representation of the equation (j3J) we consider the following: 

J_ ( y + -y _ y-y- \ = J_ (Q) 
h- { h + J yV' 1 j 

The last equation on the mesh (jSJ) can be rewritten as the mapping 

y + y(2-e 2 ) = y{y+ + y-). 

Notice that difference model @, © of equation (jHJ) can be developed from the discrete 

Lagrangian function C = y\ -qr which gives first integrals by means of a finite-difference 

analog of the Noether theorem IE] : 

h=y x 2 + ^ T = A, I 2 = + * + , X + T\ = B - ( 10 ) 

yy + 2yy + + yx{x + yx-y + ) 

Eliminating y x , x+ and y+ from the first difference integrals, we obtain the general solution 
Ay 2 = {Ax -B) 2 + l-^, (11) 

which agrees with the continuous limit up to order e 2 . 

To obtain the difference analog of the non-point transformation © we rewrite the latter as 
relations between the differentials: 

dx = y(x) 



u 2 ' u(t) 

Then, as a difference analog of the above transformation, we will use non-point transformations 
for discrete variables as follows 

1,1 1 . h + h~ 

y = -, y = — r, y =—, t = — +, r = — -, y x = -u t . (12) 

u u yy^ yy 
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where r + and r are the mesh steps in the new coordinate system. Invariant mesh (JSJ) is 
transformed into the following regular mesh in the new coordinate system: 

r + = r" = t = e, 

and the nonlinear equation @ becomes linear: 

u + — 2u + u 



TO 2 



+ u = 0. 



The general solution of the linear equation is 

u(t) = \/~Aq cos ( Bq — —t \ , 4> = arccos ( \ . (13) 



to J V 2 

The transformation (|12[) changes the first integrals (jl()j) into the following ones: 

2 , A /(f) 4>\ 1 / 2 \ 
u t + uvT = A , — t + - - - arccos — j— - cos = 5 . 

The solution family as well as the family <|13|) . are quite similar to those of the corre- 
sponding differential equations, and the superposition principle for two solutions of the difference 
scheme (jHJ), © can be formulated in the same manner as for the continuous case. 

Notice that the linearizing transformations presented above change the admitted symmetries 
of the underlining equation. The way to obtain such transformation could be viewed as follows. 
The first step is construction of an invariant difference model (that was done in 0121) with the 
conservation of Lagrangian structure, to make the Noether-type theorems applicable. Then, as 
the second step, we find a difference non-point transformation, which makes the nonlinear in- 
variant mesh regular (since otherwise even linear difference equation would be nonlinear). This 
step is quite evident from the structure of the mesh equation. The third step is to obtain a 
continuous limit of the transformation of lattice, and then the corresponding tangent transfor- 
mations in the space x, y, t, u, which preserve both constants of the two first integrals. As 
a last step, we complete the appropriate linearizing transformations for discrete models. 

It should be pointed out that the linearizing transformations are not unique. Indeed, one can 
find another relation between first integrals as far as any integral can be expressed by means 
of any smooth function of two independent first integrals. We should also add that existence of 
two explicit first integrals allows one to express the linearizing transformation explicitly from 
the space (x,y,y') to the space (t,u,u r ). In Example 4 we will deal with a generalization of 
equation (jHJ), which has two-dimensional symmetry and the only first integral. In that case we 
are able to construct linearizing transformations as connections between differentials. 

In the following examples we briefly reproduce the linearizing transformations for both in- 
variant ODEs and corresponding invariant difference equations and meshes. 



3 Example 2: ODE with 3-dimensional symmetry group 

The differential equation 

v" = Ml (14) 

has a symmetry algebra generated by the following operators 
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This equation can be obtained by the variational procedure from the Lagrangian 

which admits operators X\ and X 2 for any parameter k 7^ 1 , and therefore the Noether theorem 
yields two first integrals 

h = { -^-(y')^ +V = A°, J 2 = (k- l)(y')^ -x = B°. 

Eliminating y' we find the general solution: 

Following the same procedure as described in Example 1, one can find the non-point trans- 
formation 

x = -V(t), y = t+^(V') k + (k-l)lnV, y ' = {V') k -\ (15) 
which transforms equation (|14[) into the following linear equation 

(k - l)V" + V' = 0. (16) 
The transformation (|15|) changes the first integrals into the following ones 

(fc-l)lnV + t = A , (k - 1)V' + V = Bo. 
The latter yields the general solution of the linear equation 1)16(1 : 

V(t) = (1 - k)e^=r + B . 

As an invariant difference scheme for the equation (|14|). the following difference equations 
were constructed 0: 

2a(k-l) ( 1 . 1 n h+ h_ 



h+ + h. 



((y*) fcil - (yj feil ) = 1, 



i_ > 

1 



where the constant a should be chosen in accordance with the compatibility condition for the 
first integrals (see (2] for details). 

We again rewrite the transformation (|15|) as a connection between differentials 

and then involve transformations as a two-point change of variables 

(yx)^ 1 = -V t , (y-)^ 1 = -Vf 



X ) v ti ' 1 ) ' _ 1 

iVx)"- 1 (Vx)"- 1 



Evidently, the transformed mesh will be regular 



r+=r - =r o 



and the equation will be linear 

y _ y- y + y- 
2a(k - l)^pp- + = 0. 
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4 Example 3: ODE with a 3-dimensional symmetry group 

We consider one more three-dimensional group and its Lie algebra, 

^ d d d . , d 

*i=o-> ^2 = ^, X 3 = x— + {x + y)—. 
ox ay ox ay 

The corresponding invariant second-order ODE is the following: 

y" = exp(-y')- (17) 
It can be obtained from the Lagrangian function 
L = exp(y') + y, 

which is divergently invariant with respect to X±, X<i- The corresponding first integrals of the 
equation (fTT|) are 

exp(y') - x = A), exp(y')(l - y') + V = B . (18) 
So, the general solution of equation (|17|) is 

y = (x + B )(ln(x + B ) - 1) + A . 
The equation (|17|) can be linearized by the non-point transformation 

x = V(t), y = (l + V')lnV -V -t, y' = lnV, (19) 
which transforms the equation ()17|) into the following linear one 

V" - V' = 0. 

The transformation ()19|) changes the first integrals (|18|) into 

V'-V = A , InV' -t = B , 
which yield the general solution of the linear equation, 

V(t) = e Bo+t - A Q . 

The superposition principle for solutions of (|17jl can be established similarly to Example 1. 
The invariant difference model constructed in [2] is the following: 

A( e te_ e fe) = i ; (20) 
h+ 

hi h- , , 

— = — > 21 



where 



/3 = — j-, — - = e = const > 0, 

(l+e) 1+ - e yx 

the constants (3, e are chosen in accordance with the compatibility condition for the first inte- 
grals. 
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The above difference model preserves both all symmetries of the original differential equation 
and the Lagrangian structure (see [2]), and therefore can be integrated. The general solution is 

y = (x + B ) ln(x + B ) - (1 + e 2 )(x + B ) + A , h + = e(x + B ). 

This difference model can also be linearized. We rewrite the transformation Q19|) as 

v' = e y\ dt = ^r, 

and find the corresponding difference transformations as 

h + hr 



V f = e y * Vr = e y * 



e yx e Vx 



Applying the above transformation we linearize the difference equation 1)20(1 together with the 
mesh (j2Tj) . 

V + -2V + V~ V + - V~ . 

5 ~ = 0) T = T = T 0- 

to 2 2r 

The last fourth Lie's equation, which possesses 3-dimensional symmetry, contains an arbi- 
trary constant, and it is a first integral of some nonlinear third-order ODE (the differential 
consequence of this second-order ODE) . The linearization of this equation is rather complicated 
and is not considered here. 



5 Example 4: ODE with a 2-dimensional symmetry group 

S. Lie identified two classes of ODEs which possess a 2-dimensional symmetry group. Both 
classes contain arbitrary functions. For example, one class of equations is the following 

y" = ~F{y% 

x 

which is invariant with respect to the group with the following operators 

d d d 

ay ox oy 

It seems to be impossible to linearize the whole class of invariant ODEs, thus we consider 
just representative equations. We will choose the example which is important from the physical 
point of view. 

First, we interchange dependent and independent variables and, second, we make scaling of 
the independent variable to obtain another representation of the underlining class of equations: 

y" = y n F(x 2 y n - 1 ), 

where n is some constant, n/-l. 

Now we consider the special case of the above equation, when the function F = 1. It will be 
the general case of the power nonlinearity, i.e. the differential equation 

y" = y n , n + -1, (22) 



which possesses an isomorphic two-dimensional Lie algebra of operators 
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Notice, that this example makes sense from the physical point of view, as equation (|22|) can 
be considered as a dynamical equation for a particle in a field of the power potential. This 
equation can be obtained from the Lagrangian function 

L = y' 2 + ^-y n+ \ 
n + 1 

which admits the only operator X\. 

In contrast to Example 1, which is a special case of equation (|22|) with n = — 3, the underlining 
equation for an arbitrary n has a unique first integral: 

Ji = y' 2 - ^-y n+1 = A . 

n + 1 

The second integral can be expressed by a quadrature 
dy 



A + ^y n+1 



Bo, 

n+1 

and for the only cases 

1 , 1 - 2k 

n = - — 1, n = -, 

k l + 2k 

where k is integer, this integral can be expressed by a finite superposition of elementary functions. 
We see that equation ©, considered in Example 1, corresponds to k = — 1. In contrast to 
a general case, equation © is the only equation of power nonlinearity which possesses a 3- 
dimensional symmetry group. The exceptional potential which corresponds to equation © is 
known in theoretical physics. 

Thus, the example of the general power potential under consideration is sufficiently different 
from the above three, as we have fewer symmetries and the only first explicit integral. Never- 
theless, equation (|22j) for any n can be linearized by a non-point transformation. The only 
disadvantage is that we will express the linearizing transformation in terms of differentials. 

Again, we would like to preserve the first integral under linearization, for example, as follows: 

yl 2 -—y^=A Q = u' 2 -U 2 . 

n + 1 

Then, we can choose the change of dependent variables and derivatives 
2 



-y n+1 = u 2 , y'(x) = u'(t), 



n + 1' 

and, therefore, the transformation of differentials and derivatives could be the following 

n + 1 n i-n n + 1 1 2 

dx = | — - — | "+ 1 n™+i(it, y{x) = \ — - — |«+in"+ 1 , y(x)=u(t). (24) 

The transformation (|24j) transforms the equation ()22l) into the following linear one 

u" - u = 0, (25) 

which has the first integral 
u ' 2 _ n 2 = A . 
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Meanwhile the equation (|25|) has one more functionally independent integral 

(u' + n)e~* = .Bo, 
and therefore the general solution is 

u(t) = Ae* + Be~ t . 

Notice that one can transform backwards the above general solution and the second integral just 
for exceptional values of n. 

We will construct an invariant difference model by means of full set of difference invariants 
of the Lie algebra (|23|) on a difference stencil: 

h-i- V4- V— n-l n-1 

/ 1 = -±, h = —, h = —, h = h + y + — y—. 26 
h- y y 

As an invariant difference mesh, we will use the following relation between difference invari- 
ants: 

re — 1 n— 1 

h+y+ 4 =/i_y_ 4 , (27) 
which has the integral 

re— 1 n— 1 

h + y + 4 y 4 = e 5 e = const. 

By means of difference invariants (|26j) one can develop an invariant approximation of equa- 
tion ((221) on mesn 

( re+l re + 1 n + 1 re + 1 \ 

y + * -vj _ y> -y * )' =y+ ^. m 
h+y+ 4 h-y- 4 J n + 

As the difference analog of the transformation (|24|) . we will use the non-point transformations 
for discrete variables as follows 

71+1 1 re + 1 , .71+1.1 2+i _ ,71+1. 1 re+1 

«(*) = |— 2~ I 2 v 2 ' = I — ^ — I 2 y+ 2 , u = 1—2—1 2 v- 2 . 

I I n— 1 n — 1 n — 1 n — 1 

7" = ^ 2/ 4 ?/+ 4 , T = h y 4 y_ 4 , 

where r + and r _ are the new mesh steps. 

Invariant mesh (|27|) is transformed into the regular one 

T + = T~ = T = e, 

and the nonlinear equation (|28|) becomes linear, 

u + — 2u + u" 



TO 2 



u = 0. (29) 



Equation (|29[) has the following difference integrals 



u t 2 — im + = j4 , (n t + (3u + )q T o = Bq 



where 

_ r + yVo 2 + 4 -r + Vto 2 + 4 

^ = o ' 9 = , s ■ 

2 r + Vro 2 + 4 

The general solution of (|29|) can be found from the above integrals. Alternatively, one can 
look for a solution in the form 

u(t) = Aqi^o +Bq 2 ^. (30) 

Substituting (|3*0)l into 1)29(1 yields the needed constants 

r 2 + 2 ± 

91,2 = o ■ 
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6 Example 5: special case of an ODE 
with a 2-dimensional symmetry group 

For the special case n = — 1 the corresponding ODE 

y" = y'\ (3i) 

possesses the two-dimensional Lie algebra given by the operators 

9 ^ d d 

x i = "5"' ^2 = x— + y-^-. 
ox ox oy 

Equation (|31j) can be obtained as the Euler equation from the Lagrangian 
L = \ny+\{y'f 

which is invariant with respect to X\. 

The corresponding first integral of equation (|31j) is 

]ny-±(y') 2 = A , (32) 
while the second integral could be expressed by a quadrature 

d V =B 
V A + 2 In y °" 

Following the way of preservation of the first integral under the linearizing transformations 

\ny- l -{y'f = A, = u- l -{u')\ (33) 

we choose the change of dependent variable and derivative as follows 
\ny(x) = u(t), y'(x) = u'(t). 
Then the independent variables are connected as 
dx = e u dt. 

It follows from first integral (|33j) that the second-order ODE is the following: 

u" = 1. (34) 
Thus, equation (|31|) can be linearized by a non-point transformation 

y ( x ) = e u W, dx = e u dt, y=u, (35) 

which transforms equation (|31|) into the linear ODE (|34j) . 

The transformation (|35|) changes the first integral 1)32(1 into the following: 

u- l -(u'f = A Q . 

For the linear equation (|34l) there exists one more independent first integral 
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which yields the general solution of the linear equation in the form 

u(t) = (*±pl + A . (36) 



Thus, the linearizing transformation (|35[) yields an implicit superposition principle for equa- 
tion (|3T|) . 

An invariant difference model can be obtained from the following difference Lagrangian: 
L = ln?/ + lny + + (y, x ) 2 , 

where y x = v+ f^_ y is a right difference derivative. 

The application of the difference analog of the Noether theorem yields the appropriate 
quasi-extremal equation 

_ lny + -lny- , Q7 x 

Vx-Vx = ; — , (37) 

Vx + Vx 

where y x ~ = v ' h v is a left difference derivative, and the first integral 
ln(y+y) - y x 2 = 2A. 

The invariant difference mesh for (|37|) should be expressed by means of difference invariants 



h + = h~F 



y+ y_ h_ 

y ' y~' y 



and it is not yet defined. 

As a difference model for ()34|) we consider the following linear second-order equation 

u + — 2u + u~ 1 . , , 

— {u t -u t -) = l, (38) 



on regular lattice 

T + = r~. (39) 
The last model has the following first integrals: 

(u + + u)- u t 2 = 2A, (40) 

U *"(* + t) =jB ' (41) 

and mesh equation (|39|) has an evident integral 
r + = t~ = tq = const. 
From two first integrals (|4U|). 1)41(1 one can exclude u + and obtain a general solution 

w 2 8 

which well corresponds to the solution (|36|) of ODE (|34() . 
Now we can write down the linearizing transformation 

u = lny, u + = lny + , u^=lny~, u t = y x , uf = y x ~ , (42) 
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, lny + -lny , , lny-lny 
r = - h , r T = a , 

y + - y y-y 

which preserves both the first integral (|4U|) of equation (|38j) and the integral of a mesh. Notice 
that in the continuous limit transformation (|42|) tends to (|35j). Under the action of (|42|) the 
regular mesh ()39[) is transformed into the following non-regular invariant mesh 

In — In ~~~ 

h + —^ = h~ (43) 

y + - y y-y 

which completes the difference model for ODE (|31|) . 

Thus, invariant model (|37|). (|4"3~)) is connected with the linear model (|38|). (|3^|) by two-points 
transformation (|42|). 

The remaining class of second-order ODEs from the Lie's list contains all autonomous equa- 
tions which are invariant with respect to the only operator Xi = For this class the well- 
known procedure reduces the order of equation by one. A similar reduction can be done for any 
autonomous difference equation on a lattice which is invariant with respect to a translation of x. 



7 Concluding remarks 

It was shown that, for all examples of invariant differential equations from the Lie's list, there ex- 
ist non-point linearizing transformations (which transform the family of solutions of the original 
nonlinear ODE into the family of solutions of a linear one) and, therefore, nonlinear superposi- 
tion principle for their solutions. 

The procedure that we followed in this article is based on the preservation of at least one first 
integral for both the ODE and the difference equation, and a transformation of the invariant 
nonlinear difference mesh into a regular one. We started from the continuous case, where 
second-order ODEs possessing symmetry were already known. From we know the symmet- 
ry preserving the difference model, i.e., the difference equation and a mesh, which in addition 
conserves the Lagrangian structure of the original ODE. Then we choose the change of dependent 
variables which preserves the known first integral. These transformations can be different 
for different representations of first integrals, and we choose such change of variables, which 
makes an invariant mesh regular. Then we obtain a continuous limit of these transformations of 
a lattice and derive the corresponding transformations in the space (y,dx,u,dt), preserving all 
first integrals. For the cases of two explicit integrals (Examples 1, 2, 3), we express explicitly the 
change of variables in the space (x,y,y' ,t,u,u'). Finally, we complete the appropriate lineari- 
zing transformations for the discrete model. Notice that the developed non-point linearizing 
transformations are not contact transformation in the space (x,y,y',t,u,u'), but are such on 
the family of solutions of invariant ODEs. 

It worth mentioning that all difference and differential equations considered here are inte- 
grable (at least in quadratures) as far as they possess a symmetry algebra of the dimension not 
smaller than two. Thus, the integrability is not of interest for these equations. The result is the 
explicit non-point linearization of integrable nonlinear differential and difference models. The 
second valuable result of the linearization is a nonlinear superposition principle for solutions 
of nonlinear equations, which have definite physical meaning. Finally, we can indicate that all 
mentioned properties can be preserved in finite difference models. 
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